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ABSTRACT 

We map the spatial distribution of recent star formation over a fewxlOO Myr timescales in fifteen 
starburst dwarf galaxies using the location of young blue helium burning stars identified from optically 
resolved stellar populations in archival Hubble Space Telescope observations. By comparing the star 
formation histories from both the high surface brightness central regions and the diffuse outer regions, 
we measure the degree to which the star formation has been centrally concentrated during the galaxies' 
starbursts, using three different metrics for the spatial concentration. We find that the galaxies span 
a full range in spatial concentration, from highly centralized to broadly distributed star formation. 
Since most starbursts have historically been identified by relatively short timescale star formation 
tracers (e.g., Ha emission), there could be a strong bias towards classifying only those galaxies with 
recent, centralized star formation as starbursts, while missing starbursts that are spatially distributed. 
Subject headings: galaxies: starburst ~ galaxies: dwarf - galaxies: evolution - galaxies: individual 
(ESO 154-023, UGC 4483, UGC 9128, NGC 625, UGC 784, NGC 1569, NGC 2366, 
NGC 4068, NGC 4163, NGC 4449, NGC 5253, NGC 6456, NGC 6789, IC 4662, 
DDO 165) 



1. INTRODUCTION 

1.1. The Starburst Mode of Star Formation in Dwarf 

Galaxies 

The term starburst refers to an inte nse, finite period 
of star-formation (SF) in a galaxy (e.g. , |Gallagher||1993 



Heckman] 1997[ Kennicutt et al. 2005 1 . 'i'he increasec 
rate of ongoing Si"' in a starburst has numerous observa- 
tional signatures. The youngest, most massive stars pro- 
duce significant enhanced emission at ultraviolet (UV), 
optical, and infrared wavelengths, which can be used to 
identify regions of significant, ongoing SF. However, if 
these stars are exclusively used in the analysis, their in- 
herently short lifetimes limit identification of extended 
SF events. Indeed, a starburst can be both spatially 
and temporally larger than the present-day population of 
young, high-mass (M > 15 Mq) stars may suggest. For 
example, sites of significant, ongoing SF have been iden- 
tified Jn_jiumerous_galaxies based on their Ha emission 
(e.g. 



Lee et al. 12009 1, but the timescale for Ha emission 



5 — 10 Myr. Thus, the presence of Ha can nei- 
ther confirm nor rule out whether the SF has persisted 
on longer timescales or has been present in regions other 
than where the present-day massive stars are located. 
Moreover, SFRs measured from Ha emission are prone 
both to stochastic effects from the low number of high- 
mass stars responsible for the emission, and to increased 
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uncertainties due to the emission being a secondary pro- 
cess of the photoionization of gas in low density media 
(i.e., not from the stars themselves) ( |Lee et al.||201l| . 
SFRs measured from UV emission are found to be more 
robust, as the UV flux is emitted directly from the stel- 
lar photospheres of more numer ous, intermediate mass 
(2 M0 < M < 15 M©) B-stars ( |Lee et al.l[20lT| ). UV 
emission has the added advantage of sampling a longer 
timescale than Ha emission (~ 100 Myr versus ~ 5 — 10 
Myr) and has been show n to be an excellent tra cer of SF 



in star forming galaxies \G'\\ de Paz et al. 2007). 

Compared to UV and Ha emission, optical imaging of 
resolved stellar populations can sample higher time res- 
olutions and much longer timescales (depending on the 
photometric depth achieved by the observations). Us- 
ing high-fidelity photometry of resolved stellar popula- 
tions, the longer temporal baselines available at optical 
wavelengths can be used to reconstruct the history of SF 
over G yr-based timescales (for a review, see Gallart et al. 



2005 



and references therein) providing a much broader 
perspective of the SF experienced in a host galaxy. 

Among the many morphological types in which star- 
bursts have been seen, dwarf galaxies offer valuable lab- 
oratories in which to study the starburst phenomenon. 
Dwarf galaxies have the advantage of being both the 
most numerous morphological galaxy type (i.e., provid- 
ing a larger possible sample of local starbursts within a 
fixed volume), and having low metal content (i.e., and 
thus less internal extinction than higher metallicity sys- 
tems simplifying the analysis of the star formation his- 
tories (SF Hs)). The sohd-body rotation of dw arf galax- 
ies (e.g., Skillman et al. 1988; Skillman 1996) prevents 
the mixing of the stellar populations on short timescales, 
which facilitates spatial analysis of the rece nt SF (e.g., 
Bastian et ar]|2009[ |2011[ ). Starbursts 



Gieles et al. 2008 



in dwarf galaxies are known to create outflows of chem- 
ically enriched material into the intergalactic medium 
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( Martin et aL]|2002| |Summers et al.||2003[ [Hartwell et 
al.||2(J(J4 ). 'i'his outtlow ot gas trom a starburst has been 



linked to both potentially playing a role in transforming 
dwarf irregular galaxies into dwarf spheroida l or dwarf 
elli ptical systems (e .g., Marlowe et al. 1997| Mayer et 



al. 



2001j |van Zee et aL 20U1; Fasetto et arf|2008[ ) and 



m shaping the magnetic held of the host galaxy (e.g 



NGC 1569; [Kepley et al.||2010[ ). Finally, dwarf galaxies 
have low stellar surface densities, which makes them eas- 
ily resolved into stars, facilitating analysis of their CMDs. 

1.2. Spatially Centralized vs. Distributed Starbursts 

Historically, starbursts have been assumed to be cen- 
tralized. This assumption is based on observations show- 
ing that recent SF is often found primarily in the center 
of a starburst galaxy's optical disk. This effect is seen in 
different morp hological types, inclu ding starburst spiral 



galaxies (e.g., Caldwell et al. 1999), irregular starburst 
dwarf galaxies, and blue compact dwarf galaxies (BCDs; 
a sub-class of starburst dw arf galax ies, e.g., j Papadcros 
et al.||1996[ [van Zee et aL|[r998. ,Salzer fc Norton„1999^ . 
I'heoretical simulations oi star burst galaxies also often 
assume a central sta rburst ( e.g., Norman fc ScovilleJ1988 

99gr^ 



Mihos & Hernquist 1994; Mac Low & i-'errara 199' 



However, there is ample evidence that starbursts may 
not be a solely centralized phenomenon. For example. 



Meurer (2000) studied the UBI colors of star clusters 
and dittuse light in the starburst spiral galaxy NGC 3310 
and found that 80% of the flux in the U band was emis- 
sion from the diffuse SF throughout the optical disk of 
a galaxy and only 20% came from the centrally concen- 
trated stellar clusters. A similar study reported both cen- 
trally and non-centrally concentrated SF in ne arby BCDs 
based on the rad ial profiles oi B — V colors (Hunter fc 
Elmegreen 2006). Distributed starbursts have been re- 
ported m both low and high stellar density regions in 
individual galaxies based on resolved young stars (e.g 



Annibali et al. 2003 Cannon et al. 2003). In add! 



tion, observations have been better tit by simulations 
with widespread, shock-induced bursting SF, ra ther than 
centrally con centrated, density-dependent SF (Chien fc 
Barnes|[2010 ). The metrics used to evaluate the concen- 
tration oF'SF in these heterogeneous studies vary, pre- 
venting a direct, uniform comparison of results. However, 
the mixed findings emphasize that one cannot assume 
that starbursts are solely a centralized phenomenon. 

1.3. Re-Examining the Starburst Paradigm 

This paper is the third in a series studyi ng starbursts 
in n earby d warf galaxies. In the first paper, [McQuinn et 
al. (2010a hereafter. Paper I), we derived the star for- 
mation histories (SFHs) of twenty galaxies by matching 
the optically resolved stellar populations of the galaxies 
to synthetic color-magnitude diagrams (CMDs). These 
individual SFHs allow us to define a starburst galaxy as 
a system with unusually h igh levels of SF c ompared to 
its own historical average (Kennicutt 1998a). This per- 
spective widens the definition to include lower-intensity 
starbursts that have neither the immediate extreme ob- 
servational signatures of the youngest most massive stars 
nor the large super-star clusters / strong emission pro- 
duced by high absolute SFRs. Using this broader defi- 
nition, four galaxies in the sample were newly identified 



as recent starbursts (i.e., Antlia dwarf, ESO 154—023, 
NGC 784, and UGC 9128). Indeed, the optical images 
of these bursts do not show the larger super star clusters 
and Ha emission typically associated with well-known 
starburst galaxies such as NGC 5253, NGC 4449, or 
NGC 1569 (which are also included in our sample). Yet, 
the SFHs of these newly identified starbursts clearly show 
elevated recent star formation rates (SFRs) and qualify 
the systems as starbursts according to the criterion in 
Paper I (i.e., recent SFRs > 2 x {S¥'i\) historical) ■ Five ad- 
ditional galaxies were identified as fossil starbursts that 
have present-day SFRs decreasing to historical levels or 
below. 



In th e second paper of the series, McQuinn et al. 
(2010b hereafter. Paper II), we focused on the temporal 
behavior of the starbursts. The starburst durations were 
measured to last a minimum of 450 Myr in nineteen of 
the twenty dwarf galaxies, and as much as 1.3 Gyr in two 
of the systems. SFRs sustained on these timescales can- 
not be supported by an individual pocket of SF, due both 
to the limited availability of gas and to the likely disrup- 
tion of SF by feedback on small spatial scales. Thus, 
these longer durations suggest that a starburst can be 
more than individual pockets of flickering SF occurring 
in clusters and that a significant fraction of a galaxy may 
participate in a burst. Yet, it is unclear how much of the 
galaxy experiences a burst and whether the SF occurs in 
all regions simultaneously or propagates across a galaxy. 

In this third paper of the series, we study the spatial 
distribution of SF in fifteen of the original twenty star- 
burst dwarf galaxies with deep, wide area Hubble Space 
Telescope (HST) coverage. The goal of the paper is to 
test the assumption that starbursts are generally a cen- 
tralized phenomenon. To accomplish this goal, we trace 
the location of SF over the last few 100 Myr in each 
galaxy using positions of young, intermediate mass (2 
M© < M < 15 Mq) blue core helium burning (BHeB) 
stars from archival HST observations. In addition, we 
reconstruct the temporally resolved SFHs for regions of 
varying stellar density, thus tracking the SFR across dif- 
ferent regions of each galaxy. 

The paper is organized as follows: ^|2] summarizes the 
observations and data processing. Splpresents the metrics 
used to determine whether the starbursts are centralized 
or distributed and discusses the results of these metrics, 
fjl] summarizes the conclusions. 

2. THE GALAXY SAMPLE, STELLAR POPULATIONS, AND 
REGIONAL SFHS 

2.1. The Galaxy Sample and CMDs 

The sample of twenty starburst dwarf galaxies stud- 
ied in Papers I and II was selected from observations in 
the HST public archive. The observations were chosen 
based on three properties that ensure robust reconstruc- 
tion of the recent SFHs. First, we selected galaxies that 
lie close enough (i.e., D < 5 Mpc) to resolve their stel- 
lar populations with HST's imaging i nstruments — the 



Advanced Camera for Surveys (ACS; Ford et al 



or the Wide Fiel d Planetary Camera 2 (WFPC2; |Holtz 



1998^ 



man et al. 1995). Second, we required the archival ob 
servations to include both V and I band images of a 
galaxy. Third, we selected systems where the I band 
observations' photometric depth reached a minimum of 
~ 2 mag below the tip of the red giant branch — a re- 
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quirement for accurately constraining the recent SFH o f 



NGC 784 



2002 Dohm-Palmer et al. 2002) 



the galaxy ( |Dolphin| ^ 

One galaxy (SBS 1415+437) did not meet this final re- 
quirement but was included in the original studies for 
comparison purposes. It is not analyzed in this paper. 
While our archival-defined sample spans a wide param- 
eter space for starburst dwarf galaxies including many 
of both the stronger, well-known starburst dwarf galax- 
ies in the nearby universe (e.g., NGC 1569, NGC 4449, 
NGC 5253) and the smaller, less well-known starbursts 
(e.g., NGC 4068, UGC 9128), it is not a volume limited 
sample. 

From this group of nineteen galaxies, we applied two 
additional cuts to ensure global spatial analysis was pos- 
sible. First, we required that the HST observations cov- 
ered at least 50% of the deprojected area out to the 
B-band 25 mag arcsec"^ isophotal di ameter (D25), as 
defined by [Karachentsev et al. ( 2004 ) . This criterion 



eliminated three galaxies with limited areal HST cover- 
age (NGC 4214, NGC 6822, Holmberg II). Second, we 
required that the number of core helium burning stars 
allowed for a meaningful spatial analysis, which required 
that the amplitude of the burst must be greater than 
10"'^ Mq yr"-'^. This criterion eliminated one galaxy 
(Antlia) from the sample. The final sample consists of 
the fifteen dwarf galaxies listed in Table [l] Also included 
in Table [1] is the percentage of the deprojected D25 area 
covered by HST and the equivalent area in kpc^ . The de- 
projected area was calculated by correcting the angular 
area of each galaxy for the optical inclination angle (i), 
Galactic extinction in the B-band, and the distance to 
the system. The specific HST instruments, observation 
times, and H ST program I Ds ca n be found in Paper I 
(Table 1), and | Weisz et~al] ( [2008| . 

All photometry wa s measureTby DOLPHOT or HST- 
phot ( [Dolphin 2000) from the standard HST pipeline 
processed and cleaned images. The data sets were pro- 
cessed uniformly, ensuring reliable comparison between 
galaxies. Artificial star tests were performed to mea- 
sure completeness using the same photometric software 
packages listed above. The details of the photometric 
processing for fourteen of these galaxies are presented 



in Paper I and for one galaxy (DDO 165) in Weisz et 
al. (20081. CMDs for all galaxies can be found in I' lg- 
ure 2, Paper I. Here in Figure [l] we present an example 
CMD for NGC 784, whose data reach our minimum re- 
quired photometric depth. Representative uncertainties 
are shown per magnitude by crosses in the left hand side 
of the CMD. 

Stellar evolutionary features that are readily apparent 
in the example CMD of Figure [T] include the main se- 
quence (MS), blue and red core helium burning (BHeB 
and RHeB respectively) sequences, the red giant branch 
(RGB), and the asymptotic giant branch (AGB). The 
upper MS, the BHeB, and the RHeB stars are populated 
by stars > 2 Mq and thus ar e unambiguous indicators 



of recent star form ation (e.g., Gallart et al. 

and references therein). 



2005 Mc- 



Quinn et al. 2011 and references therein). Stars with 
a mass of 2 Mq or greater will evolve off the MS and 
begin burning helium in their cores. These RHeB stars 
form a sequence slightly hotter than the RGB stars. At 
brighter magnitudes, the RHeB stars are well-isolated 
from other populations, but at fainter magnitudes they 
begin to become indistinguishable from the RGB stars. 
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Figure 1. The CMD of NGC 784 showing the typical photometric 
quality of the data used in the analysis. Photometric uncertain- 
tics arc shown in the left hand side of the CMDs. Major stellar 
evolutionary stages are marked on the CMD. Contours are used 
to show the photometric structure of the stellar populations in the 
red giant branch and lower main sequence. 



The RHeB stars remain at cooler effective temperatures 
until their opacity changes due to internal convective 
mixing. These stars then become hotter and rapidly tra- 
verse the "blue loop" in the CMD becoming BHeB stars, 
locat ed just redward of the up per MS stars in CMD space 



(e.g., Stothers & Chin 1991). When the core helium is 



exhausted, the stars will evolve further into either AGB 
stars, moving redward in the CMD to cooler tempera- 
tures, or, erupt as supernovae. 

2.2. HeB Stars Trace the Locations of Recent SF 

As described in §2.1| and seen in Figure [T] the BHeB 
and RHcB stars occupy a unique space in the CMD. 
Along the RHeB and BHeB sequences, there is a nearly 
one-to-one correspondence between the luminosity of a 
HeB star and the mass of the star. The masses and 
stellar ages of each point on the sequences can be es - 
timated using stellar evolution isochrones (see j |2.2.2 ). 
Thus, HeB stars in a given luminosity range can be ro- 
bustly age-dated using st ellar evolution isochrones (e.g., 
Dohm-Palmer et al. 19971, and the spatial distribution of 
HeB stars of different ages can be used to trace the sites 
of recent SF. Note that the luminosity of a HeB star does 
undergo small changes as the star evolves. This minor 
change in luminosity is smaller than our photometric un- 
certainties and we do not correct for it. Here we present 
our methodology for isolating HeB stars, dating them us- 
ing stellar evolution libraries, and thus tracing the sites 
of recent SF. We use one galaxy (NGC 784) as an ex- 
ample. In the Appendix, we present the results of this 
analysis for the entire galaxy sample. 

2.2.1. Isolating the HeB Stars in an Optical CMD 

At low metallicities, both blue and red HeB stars are 
relatively isolated in an optical CMD as seen in Figure [T] 
The bright BHeB stars form a sequence that is approxi- 
mately parallel to and slightly redward of the upper MS; 
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the exact shape of the BHeB sequence depends on the 
host galaxy's metaUicity, and to a lesser extent, the SFH. 
At fainter magnitudes, the sequence of BHeB stars trends 
redder and eventually blends into the red clump. The 
typical intrinsic V— I color separation between the MS 
and BHeB stars is 0.1 — 0.3 mag. This distinction can 
become blurred in galaxies with higher metallicities, in- 
ternal reddening, and photometric crowding. However, 
the starburst galaxies in our sample have low metaUicity 
and low internal reddening. They also have less stellar 
photometric crowding than larger spiral galaxies on aver- 
age, although dwarf starburst galaxies do have regions of 
higher stellar density. Overall, photometric uncertainties 
and crowding are relatively small and do not significantly 
blur the MS and BHeB populations, as can be seen in the 
representative CMD in Figure[T] The separation between 
the populations is more apparent at higher luminosities, 
and becomes less distinct at lower luminosities depend- 
ing upon the photometric uncertainties. There are three 
exceptions in our sample (NGC 1569, NGC 4449 and 
NGC 5253), where the MS and BHeB populations begin 
to blend due to photometric crowding and differential 
extinction. Nonetheless, we include these three galaxies 
in the analysis noting higher uncertainties in the spatial 
analysis of BHeB stellar population. Note also the RHeB 
stars form a sequence blue- ward of the RGB stars; these 
stars blend with the RGB stars at luminosities below the 
tip of the red giant branch (TRGB) making them diffi- 
cult to uniquely separa te in a CMD for ages > 100 Myr 



for most galaxies ( ,McQuinn et al.| 2011) 
The technique of separatin g the VlS, BHeB, and RHeB 



stars is de scribed fully in Dohm-Palmer fc Skillman 



|2002) and [McQuinn et al.|p01ip . While the MS and 



eB sequences are readily identifiable in the CMDs, 
separating these sequences by eye can be ambiguous, par- 
ticularly at lower luminosities (see Figure [l] as an exam- 
ple). To robustly separate the MS and BHeB stars, we 
created histograms of the V— I colors of all stars in the 
MS, BHeB and RHeB regions of the CMD, in intervals 
of 0.5 to 1.0 magnitudes in My, depending on the pho- 
tometric uncertainties and on the number of HeB stars 
present. Binned in this way, the concentration of stars 
forming the MS, BHeB, and RHeB sequences form three 
distinct peaks in the stellar color distribution. These 
three histogram peaks mark the center position in the 
CMD of each stellar sequence. We conservatively se- 
lected the stars with colors of ±0.1 mag of the peaks for 
both the blue and red HeB populations in each galaxy. 

At fainter magnitudes, the BHeB and RHeB stars be- 
gin to merge with the red clump and RGB stars respec- 
tively, creating a lower luminosity limit for selecting HeB 
stars from a CMD. This limit is specific to each data set, 
because it depends on the metaUicity, photometric uncer- 
tainties, and SFH of each galaxy. The individual lower 
limits were determined by creating synthetic CMDs from 
the SFHs (Paper I) and by mimicking photometric com- 
pleteness and observational uncertainties using artificial 
star tests performed on the original data sets. In the 
left panel of Figure [2] we plot age versus luminosity for 
the stars in the BHeB region of a synthetic CMD, based 
on the best fit SFH of NGC 784. The one-to-one age- 
luminosity relation of the BHeB stars is apparent in Fig- 
urej2] until an age of ~ 200 Myr, where there are stars of 
multiple ages present for a given luminosity. It is at this 
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Figure 2. The age distribution of stars in the bhie (left panel) 
and red (right panel) HeB stellar regions as a function of magnitude 
for a simulated CMD matching the completeness of the data for 
NGC 784. There is a one-to-one correspondence between age and 
luminosity for the HeB stars until the HeB stars begin to mix with 
red clump or RGB stars at fainter magnitudes. This luminosity is a 
lower limit for unambiguously selecting HeB stars from other stellar 
populations in a CME). For this galaxy, the limit corresponds to 
an age of ~ 200 Myr for BHeB stars and ~ 150 Myr for the RHeB 
stars. The limit is unique for each system and dependent upon 
photometric uncertainties, crowding, metaUicity, and the SFH of 
the host galaxy. 



magnitude where the BHeB sequence for the simulated 
CMD of NGC 784 becomes mixed with older red clump 
stars. BHeB stars cannot be unambiguously selected at 
magnitudes fainter than this limit. Likewise, as seen in 
the right panel of Figure [2l the RHeB region becomes 
contaminated with older RGB stars at low luminosities, 
corresponding to an age limit of ~ 150 Myr. In general, 
we find that the BHeB stars can be traced further back in 
time than the RHeB stars, making them a longer range 
chronometer of recent SF. 

2.2.2. Dating the HeB Stars using Theoretical Isochrones 

Once the blue and red HeB stars were identified in the 
CMD, we used stellar evolution isochrones to age date 
the stars. Stellar evolution isochrones track the paths of 
stars of differ ent mass in a CMD. We use the theoretical 
isochrones of Girardi et al. (12000) for stars with masses 
< 7 M0, with improved treatment of the asympt otic gi- 



ant branch (AGB) stars added by Marigo et al. (|2008 



At higher ma s ses (M >7 M©), we use the isochrones of 



Bertelli et al. (1994). This set of isochrones, collectively 
known as the Padova stellar evolution isochrones, is used 
to age-date the HeB in the observed CMDs. Isochrones 
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Figure 3. The CMD for NGC 784 overlaid with stellar evolution 
isochrones from 40 Myr to 160 Myr in increments of log St = 
0.2. The blue and red HeB stages of evolution, easily identified 
in the isochrones, are in general agreement with the corresponding 
sequences of HeB stars in the observed CMD. 



were matched in metallicity and adjusted for both the 
internal and foreground extinction for each galaxy. The 
individual metallicities were taken from H ii region abun- 
dance measurements in the literature whenever possi- 
ble. In five cases, the metallicity was estimated from 
the luminosity-metallicity relation for dwarf g alaxies (i.e, 
Zaritsky et"ar]p94l [Tremonti erar][2004l |Lee et al. 
2006p . The Ay extinction measurements were estimated 
from the SFHs (Paper I) and transforme d to Aj extinc - 
tion using the York Exti nction Solve r tool fMcCall 2004 ), 
assuming the Cardelli et al. ( 1989 1 reddening law with 
R= 3.1. These values are m good agreement wit h the 
values from the dust maps of Schlegel et al. (1998), but 
also incorporate estimates of internal extinction. The 
ad opted metallicity and extinction values can be found 
in|McQuinn et al.|(!20n| Table 1). We opted not to use 
the isochrones incorpo rating circumste llar dust models 



p^ressan et al.^l998;,Groenewegen|2006 ), as they primar- 
iiy attect the AUB stage ot evolution and do not affect 



the location of HeB stars in the CMDs. 

In Figure [31 we present a representative comparison 
between the the stellar populations for NGC 784 and the 
stellar evolution isochrones up to 160 Myr. The blue and 
red HeB stages of evolution are easily identified in the 
isochrones, as stars in these stages occupy extreme blue 
and red turning points in the isochrones. The isochrones 
are in general agreement with the MS locus and the blue 
and red extremes of the blue loop sequence. 
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Figure 4. The luminosity functions of the BIfeB (top panel) and 
RHeB (bottom panel) stars for NGC 784. Listed above the his- 
togram functions are the ages of the stars derived from Padova 
stellar evolutionary models. The luminosity functions are plotted 
only to the fainter magnitudes below which the HeB stars begin to 
merge with other stellar populations (i.e., Mv ~ —3.5 for the BHeB 
stars and ~ -3 for the RHeB stars). The BHeB and RHeB 
luminosity functions are similar in profile, as expected. Their spa- 
tial distributions, shown in Figure [5] are also very similar. Thus, 
we focus on only the BHeB stars for the remainder of the analysis. 

In Figure |4j we plot the luminosity function of the 
BHeB and RHeB stars selected from the CMD for 
NGC 784, along with the ages derived from the stellar 
evolutionary isochrones. As expected, there are smaller 
numbers of short-lived high mass stars. Uncertainties in 
the luminosity function include both Poisson errors and 
incompleteness measurements. The luminosity functions 
extend to just above the lower luminosity limit where the 
HeB sequences begin to merge with other stellar popula- 
tions in the CMD. The luminosity function of the RHeB 
stars is similar to that of the BHeB stars providing a 
consistency check on our technique. 

2.2.3. Mapping the Distribution of HeB Stars at Different 

Ages 

We tag each bin in the luminosity function with an 
equivalent age, determined from the stellar isochrones. 
Representative ages are shown as arrows in Figure|4] The 
locations of the HeB stars in each age bin can then be 
mapped directly onto the optical image of the galaxy. We 
chose to map the HeB stars at four different luminosity 
(or age) bins. The bin selection was different for each 
system and was based on the structure of each observed 
CMD, the number of HeB stars present, and the lumi- 
nosity at which t he Hc B sta rs me rged into other stellar 
populations (see ^2.2.1 and ^2.2.2|. 

In Figure [5l we plot the positions of the BHeB and 
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Location of BHeB Stars 



Location of RiHeB Stars 











Figure 5. The location of the BHeB stars (left panels) and RHeB 
stars (right panels) of different ages in NGC 784 overlaid on a 
smoothed V band llST image. The areas of the asterisks are scaled 
by the average SFR in that age bin to account for the IMF and 
stellar lifetimes so that the total area in the asterisks is proportional 
to the SFR in that age bin. Thus, while the number of points in a 
given age bin is determined by the number of stars, the total area 
of the points represents the SFR after normalization by the IMF. 
This allows the simultaneous presentation of the location of the 
SF, the total SFR, and the number of stars from which the SFR is 
computed. Similarly to the luminosity functions in Figure |4] the 
distribution and number of RHeB stars is in good agreement with 
the BHeB stars. We present only the distribution of the BHeB 
stellar population for the remainder of the paper as they can be 
uniquely separated in a CMD farther back in time. 



RHeB stars on the optical image of NGC 784 for each 
of four age bins. The background image in each panel is 
a smoothed V band image made by applying a 25 pixel 
(i.e., 1.25") radius Gaussian kernel to the original HST 
observation. Included in Figures |4] and [5] are the results 
for both the blue and red HeB stars in NGC 784. Both 
the luminosity functions and spatial patterns in the two 
stellar populations are similar, as expected. Thus, we 
focus on the BHeB population for the rest of the paper 
as these stars are both more numerous than the RHeB 
stars (e 



Dohm-Palmer fc Skillman||2002[ McQuinn et 
al.|2011[ )" a^ can be traced to fainter magnitudes in the 



observed CMDs (i.e., older ages). 

For each age bin, the HeB stars are plotted as white 
asterisks in Figure [5) The areas of the asterisks are 
scaled by the average SFR in that age bin to account 
for the IMF and stellar lifetimes so that the total area 
in the asterisks are proportional to the SFR in that age 
bin. Thus, while the number of points in a given age 
bin is determined by the number of stars, the total area 



of the points represents the SFR after normalization by 
the IMF. This allows the simultaneous presentation of 
the location of the SF, the total SFR, and the number 
of stars from which the SFR is computed. As there will 
naturally be fewer high mass, younger stars, scaling the 
asterisks in this way enables a more direct comparison 
of the amount of SF occurring in each age bin. For ex- 
ample, in Figure |5j the total area in the asterisks in the 
most recent — 25 Myr compared to those 25 — 50 Myr 
ago indicates that the SFR was higher more recently. 

Combining the above analysis, in Figure |6j we plot the 
CMD, BHcB luminosity function, and maps of the spa- 
tial distributions of the BHeB stars in four age bins for 
NGC 784. Similar figures mapping the sites of recent SF 
in each galaxy of the sample can be found in the Ap- 
pendix. Two of the galaxies, NGC 6789 and NGC 4163, 
host fossil bursts than ended 80 Myr and 100 Myr ago 
respectively. Similarly to the other galaxies in our study, 
we plot the spatial distributions of the BHeB stars in 
these galaxies in four time bins, but define the first time 
bin in the maps to be the post-burst time period. We 
limit our subsequent spatial analysis in these two galax- 
ies to the SF occurring during the active starburst (i.e., 
prior to this post-burst time period). 

2.2.4. Inherent Age Limit and Diffusion Timescale of Stars 

Using BHeB stars to trace the sites of recent SF is 
based on the assumption that the stars have not migrated 
significantly from their birth sites. Dwarf galaxies are 
known to be solid-b ody rotators lacking sheer and differ - 
ential rotation (e.g., Skillman et al. 1988[ Skillman|1996 ). 
Their dynamics implies that there should be little radial 
mixing of stars over timescales considered here. The lim- 
iting timescale for our technique is therefore the diffusion 
timescales for dissolving SF structures. 

Stellar populations within galaxies have been shown to 
be born with a high degree of substructure made up of 
a conti nuous distribution of clusters, groups, and associ- 
ations ( [Bastian et al. 2010 and references therein ) . For 
a sample of dwart galaxies, Bastian et al. (2011) found 
these structures persist on the order ot hundreds of Myr. 
This timescale was measured down to the limiting pho- 
tometric depth of the data and is thus a lower limit. 
These findings are similar to previously reported lower 
limit lifetimes of ~ 80 M yr for stellar structures in the 
SMC (Giclcs et al "20081 and - 175 Myr for the LMC 
( [Bastian et al. 2009| ). While the stellar structures could 
not be probed on longer timescales given their data, it 
is assumed that the stellar structures will continue to 
dissolve, including some radial migration of the stellar 
populations. Although our data are frequently sensitive 
to HeB stars as old as 400 Myr, we limit the spatial 
analysis to HeB stars younger than 250 My r , the ap- 



proximate lower limit found by Bastian et al. (20111 for 



stellar structure to be maintained in dwart galaxies. We 
list the adopted limiting timescale for each galaxy in Ta- 
ble [H 

3. EVALUATING THE SPATIAL CONCENTRATION OF 
THE STARBURSTS 

As described in |T] while starbursts have often been 
assumed to be centrally concentrated, there is also ob- 
servational evidence that starbursts can be distributed in 
nature. The current data set of temporally and spatially 
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Maps of BHeB stars of varying ages on smootfied HST images 

Figure 6. Top left panel: CMD of NGC 784 with the region con- 
taining BHeB stars used in the analysis marked in blue. Top right 
panel: The luminosity function of the BHeB stars identified in the 
CMD. Middle and lower panels: The location of the BHeB stars of 
four different age ranges overlaid on smoothed V band HST images 
of the galaxy. The area of the points is scaled by the SFR enabling 
a direct comparison between age bins. The white box drawn on the 
first image encompasses the area used in reconstructing the SFH 
of the central region. The region outside of this box comprises 
the outer region used in our analysis. Note that the HST images 
are not shown in WCS orientations. Similar figures for the entire 
sample can be found in the Appendix. 

resolved stellar populations allow us to measure the spa- 
tial extent of the SF. We focus on quantifying this spatial 
distribution using two timescales. The first is over the 
last 100 — 250 Myr using BHeB stars, and the second is 
over the duration of each starburst event. In the follow- 
ing subsections, we discuss three techniques employed to 
evaluate the concentration of bursting SF over these two 
timescales. 

3.1. Three Metrics for Measuring the Concentration of 

SF 

We employ three techniques to measure the spatial 
concentration of the SF. The first technique compares 
the spatial distribution of the young stellar populations 
as traced by the BHeB stars (ages < 100 - 250 Myr) 
with the underlying older stellar populations as traced 
by the RGB stars (ages > 1 Gyr). This metric evaluates 
the concentration of SF over one timescale corresponding 
to the oldest BHeB stars (see Table [ij Column 6). The 
second technique compares the SFHs reconstructed from 
the central and outer regions. The third technique esti- 
mates the fraction of stellar mass created in the central 



region. Because both the second and the third metrics 
use the temporally resolved SFHs, they can be evaluated 
over different timescales. We chose two timescales. The 
first matches that defined by the BHeB stars used in the 
first metric. The second is equal to the duration of the 
burst events (see Table [l] Column 5). 

3.1.1. First Metric: BHeB Star Concentration 

Our first technique uses the BHeB stars selected in 
pXl| shown in Figure [6] for NGC 784 and in the Ap- 
pendix for the entire sample. We compare the distri- 
butions of BHeB stars to those of the RGB stars. The 
RGB stars were selected from the well-populated and 
well-defined RGB based on magnitude and color cuts. 
For example, in Figure |6j we selected RGB stars within 
a region defined by -3.2 <My < -1 and 0.8 <V-I< 1.3. 
The actual region used is curved and follows the distribu- 
tion of RGB stars in the CMD. We fit the distributions 
of the BHeB and RGB stellar number densities with 2-D 
elliptical Gaussian functions. While the stellar distribu- 
tions are not perfectly described by a Gaussian profile, a 
Gaussian function provides a good first order approxima- 
tion of the distribution. Uncertainties are estimated to 
be of order 10% in fitting the full-width at half-maximum 
(FWHM) of the stellar distributions based on the er- 
rors on the fit parameters. One system, NGC 6789, was 
dropped from this particular analysis due to difficulties 
fitting low numbers of BHeB stars. The degree of central 
concentration of the young stars can be measured using 
the ratio of the FWHM of the two stellar populations' 
distributions. 

In Figure [Tj we plot the cumulative distribution of this 
ratio for the sample. There is a continuous spectrum of 
the spatial extent of recent SF in these systems. There is 
no dependence on the depth of BHeB photometry in the 
distribution of the ratio measurements. In some galax- 
ies, the BHeB stars are confined to within 25% of the 
RGB spatial distribution, indicating centralized SF. In 
other galaxies, the BHeB stars extend to more than 80% 
of the RGB spatial distribution, indicating distributed 
SF. However, there are a number of galaxies where the 
BHeB stars extend over roughly half the RGB spatial 
distribution, emphasizing the point that the spatial dis- 
tribution of SF lies along a continuum. The SF in these 
galaxies are neither clearly centralized nor clearly dis- 
tributed, but fall in an intermediate category. We will 
discuss these results further when comparing the results 
of the three metrics. 

3.1.2. Second Metric: b Parameter ratios 

Our second technique for measuring the concentration 
of SF analyzes the SFHs of the central and outer re- 
gions of each galaxy. We define the "central" regions 
to be within where the V band surface brightness drops 
to 10 — 25% of the peak flux in the unsmoothed HST 
images. In one case, NGC 2366, a large star forming 
complex lies outside the central region making the flux 
cut-off between central and outer fall at the upper end of 
this range. In Figure |6j the central region for our exam- 
ple galaxy NGC 784 is outlined on the smoothed images 
of the galaxies; the outer regions encompass the remain- 
der of the fleld of view. We used rectangular regions to 
deflne the central regions. We found that our analysis is 
robust against small changes in the delineation between 
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Ratio of the Spatial Extent of BHeB to RGB Stars 

Figure 7. Cumulative distribution of the ratio of the spatial ex- 
tent of the young BHeB stars (< 100 - 250 Myr) versus older RGB 
stars (> 1 Gyr). The distribution shows a continuous spectrum of 
the spatial extent of recent SF in the sample. In some of the galax- 
ies, the FWHM of the BHeB stars spatial distribution is confined 
to within 25% of the FWHM of the RGB distribution, indicating 
highly centralized SF. In other systems, the FWHM of the BHeB 
stars extends to over 80% of the RGB distribution, indicating dis- 
tributed SF. There are also a number of galaxies where the BHeB 
stars extend to roughly half the extent of the RGB stars, empha- 
sizing that the spatial distribution lies along a continuum. 



"central" and "outer" regions and the shape of the re- 
gions, because of the large number of stars. The central 
and outer regions for the entire galaxy sample are shown 
in the Appendix. 

We reconstructed separate SFHs for the central and 
outer regions, using the same technique and method- 
ology that we applied to the whole field of view in 
Paper I. Briefly, we derived the star format ion rate 
(SFR jt, Z)) applying a CMD fitting technique CDolphinI 
2002 1 w ith stellar evolutionary models from Marigo & 
(jirardi| ( 2007 ) to the photometry and artificial star lists 
in the central and outer regions. When necessary, addi- 
tional artificial stars were generated to adequately sam- 
ple the completeness in small areas. We fit for distance, 
extinction, and stellar age, and constrain the metallic- 
ity to increase with time (except in NGC 2366 where 
the photometry reached a full magnitude below the red 
clump). The fit distances were in good agreement with 
TRGB measured distances from the literature. We refer 
the reader to Paper I for a full discussion of the SFHs, 
the time binning used, and uncertainties. 

In Figure [Sj we present the SFHs of the last 1.5 Gyr 
for both the central and outer regions in the galaxy sam- 
ple. The CMD fitting technique used to reconstruct these 
SFHs uses information from all the stellar populations 
in the CMDs, not just the HeB stars discussed in the 
preceding analysis. Thus, the SFHs are not limited by 
the same constraints in selecting blue or red HeB stars 
and can extend over a much longer timeline. A horizon- 
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Figure 8. The SFHs of regions of the central and outer regions 
of the galaxies over the past 1500 Myr. The central regions are 
defined for each galaxy in Figurc [T4] in the Appendix. The duration 
of the starbursts from (Paper II) are drawn as a bar at the top of 
each plot. The lifetime average SFR (i.e., birthrate parameter= 1) 
derived from each region is drawn as a solid gray line. 



tal bar shows the burst durations for the each galaxy as 
measured in Paper I. 

A useful metric to quantify starbursts is the birthrate 
parameter, defined as the ratio of the current SFR to 
the past aver age (b 

Starbursts have been 



nicutt 



1998a). 



SFR / (SFR), [Scalol [1986 
found to 



Ken- 
have 

birthrate parameters greater than 2; for comparison, cur- 
rent birthrate parameters for normal star-form ing galax- 
ies are typically of order 0.5 ( Kennicutt et al.„ 2005). In 
Figure [U the solid horizontal Ime designates a value b= 1 
for the individual region. 

We are interested in assessing whether the SF in the 
central regions has been more enhanced than the SF in 
the outer regions. Thus, wc compare the ratio of the 
birthrate parameter in the central and outer regions for 
the sample. In Figure |9j we plot the cumulative distribu- 
tion of this ratio, for the two timescales we selected. In 
the most recent 100 — 250 Myr (solid blue line), some of 
the galaxies have ratios greater than 10, indicating the 
central regions showed more enhanced SF than the outer 
regions (i.e., highly centralized SF). However, nearly half 
of the galaxies have birthrate parameters in their outer 
regions within a factor of two of the central regions. This 
indicates that both the central and the outer regions ex- 
perienced enhanced SF activity. This reinforces our re- 
sult from the first metric that shows a continuous spec- 
trum of the spatial extent of recent SF. On the longer 
timescale of the burst durations (dashed red line), the 
distribution in the ratios is much smaller. This tighter 
range in the ratio of the birthrate parameter is partly due 
to the differences in the time binning used in the most 
recent 100 — 250 Myr and that used over the longer du- 
ration timescales, as averaging the birthrate parameter 



Spatial Distribution of Star Formation 



9 



30 - 
20 r 
10 - 



30 - 
20 

10 - 



40 - 
20 - 



60 - 
40 - 
20 - 



NGC 6789 Central 



NGC 6789 Outer 



NGC 4068 Central 



NGC 4068 Outer 



IC 4662 Central 



IC 4662 Outer 



1 DDO 165 Central 




1 DDO 165 Outer 



1500 1000 500 
Age (Myr) 



Figure [sj Cont. 




1000 500 
Age (Myr) 



60 M 
40 ^ 
20 i 



NGC 5253 Central 



400 - 
200 - 



NGC 5253 Outer 



NGC 1569 Central 



g 200 
100 



300 r 
200 - 
100 - 




NGC 1569 Outer 




- 1500 1000 500 

Age (Myr) 



1500 1000 500 
Age (Myr) 



Figure [s] Cont. 



200 
150 
100 
50 



100 
50 



" 200 
E 150 



; ESO 154-023 C 


entral 


Duration ; 
ll**4s!S(i| 


: 1 1 1 1 1 1 1 1- 

; ESO 154-023 C 


III. 
)uter 


: 1 1 1 1 1 1 1 


mm: 


r NGC 2366 Cent 


ral 


[) 




_ 1 1 1 1 1 1 1 |T 
r NGC 2366 Oute 

F- r , ' , ^ 


r 



.1111,111 

L NGC 625 Centr 


III. 
al J 


DiirrUion . 


. 1 r — ] — 1 — M — ] — ' 1 ' 1 r 1 ' 1 . 
L NGC 625 Outer J 


: , , , , 1 , , , , ^ 


r NGC 784 Centr 


al H 




; 1 1 1 1 1 1 1 1 h 

r NGC 784 Outer 




- . . . . 1 . . . 





100 

80 

60 

40 

20 

100 

80 

60 

zn 

40 3 
20 ~ 



1000 500 
Age (Myr) 



1000 500 

Age (Myr) 



Figure [s) Cont. 



over longer timescales eliminates peak values. However, 
despite this bias, there is still a significant range in the 
spatial distribution of SF in these galaxies when mea- 
sured over the duration of the starbursts. 

3.1.3. Third Metric: Fraction of Stellar Mass Formed 

Our third technique estimates the fraction of the stel- 
lar mass formed in the central region. The stellar mass 
was estimated by integrating the SFRs in the central and 
outer regions over the same two timescales used in com- 



paring the central and outer SFHs. Note that this third 
metric is not strictly independent of the second metric as 
they both use the SFHs to measure the spatial concentra- 
tion of SF. However, the ratio of the birthrate parameter 
in the second metric compares how intense the SF is in 
the inner and outer regions relative to their historical av- 
erage, whereas the third metric integrates the SFH only 
over a specific time period. In other words, these two 
metrics both use the SFHs but have different normal- 
izations and thus provide different measurements of the 
spatial conc entr ation of SF. 

In Figure |10[ we present the cumulative distribution 
of the fraction of stellar mass formed in the central re- 
gions. Similarly to the other two metrics, in the most 
recent 100 — 250 Myr (solid blue line) the fraction of 
stellar mass formed in the central regions shows a con- 
tinuum of values. In some galaxies, more than 80% of 
the stellar mass was formed in the central regions, indi- 
cating highly centralized SF. In other systems, less than 
25% of the stellar mass was formed in the central re- 
gions, indicating distributed SF. However, there are also 
a number of galaxies where roughly half of the stellar 
mass was formed in the central regions, emphasizing the 
point that the spatial distribution of SF lies along a con- 
tinuum. When measured over the durations of the bursts 
(dashed red line), the distribution of the ratio of central 
to outer stellar mass in Figure 10 is narrower. No galaxy 
is forming more than 80% of their mass in the central 
region, indicating that the SF is less centralized over this 
longer timescale. 

3.2. Are the Bursts Centralized or Distributed: A 
Comparison of Three Metrics 

In Figure [TT| we compare the three metrics for evaluat- 
ing the concentration of SF over the most recent 100—250 
Myr. We mark the highly centralized SF end of the con- 
tinuous spectrum with a red "C", and the distributed 



10 



McQuinn et al. 



Distributed 



Centralized 



Over the Duration of Burst 

J 



Most Recent 100-250 Myr 



5 10 15 

Ratio of Central to Outer b Values 

Figure 9. Cumulative distribution of the ratios of birthrate 
parameters in the central to outer regions evaluated over two 
timescales. The ratio at more recent times (solid blue line) shows 
a larger range in the distribution. While a number of galaxies have 
high ratio values indicating highly centralized SF, nearly half have 
birthrate parameters in the outer regions comparable to within a 
factor of two of the values in the central regions, indicating greater 
spatial distribution in the SF. Not all starbursts are highly central- 
ized. Over the longer timescale of the burst durations (dashed red 
line), the ratios are smaller overall. This smaller range is partly 
due to the differences in the time binning used in the most recent 
100 — 250 Myr and that used over the longer duration timescales. 
However, despite this bias, there is still a significant range in the 
spatial distribution of SF in these galaxies when measured over the 
duration of the starbursts. 



SF end with a red "D" . The results for individual galax- 
ies in each metric are listed in Table [2j Two galaxies 
(NGC 4163 and NGC 5253) had no detectable SF oc- 
curring in the outer regions over this timescale and are 
not plotted. Based on the SF occurring in the central 
regions in these two galaxies, the ratio of the BHeB to 
RGB stars spatial extent were 21% and 14%, indicating 
highly centralized SF. 

As seen in Figure [TT] the three metrics are corre- 
lated. The systems with much smaller distributions of 
the younger BHeB stars relative to the older RGB stars 
show lower birthrate parameters in the outer regions of 
the galaxies relative to the central regions. These sys- 
tems also have a larger fraction of stellar mass being 
formed in the central, high surface brightness regions. 
Likewise, the galaxies where the spatial distributions of 
younger BHeB stars are roughly equivalent to the spatial 
extent of the older RGB stars also show birthrate param- 
eters in the outer regions that are roughly comparable to 
the central region values. Additionally, smaller fractions 
of the stellar mass are formed in the central regions in 
these galaxies. 

As seen in Figure |11[ there is a continuum in the spa- 
tial extent of the starbursts. Some galaxies, such as 
NGC 4163, NGC 5253, and NGC 6789, do show highly 
centralized SF in the most recent 100 — 250 Myr. In con- 



O 



3 



Distributed 



Over the Duration of Burst - 



Centralized 



Most Recent 100-250 Myr 



0.2 0.4 0.6 O.B 1 

Fraction of Stellar Mass Formed in Central Region 

Figure 10. Cumulative distribution of the fraction of stellar mass 
formed in the central regions of the starbursts over the most recent 
100 — 250 Myr (solid blue line) and over the lifetime of the burst 
(dashed red line). Similar to the other two metrics, a continuous 
distribution is seen in the figure, particularly at more recent times. 
In some galaxies, over 80% of the stellar mass is formed in the 
central regions, indicating centralized SF. In two systems, only 
25% of the stellar mass is formed in the central regions, indicating 
distributed SF. Over the duration of the bursts, there are more 
galaxies with a greater percentage of their stellar mass formed in 
the outer regions. 



trast, other galaxies, such as NGC 4068, ESQ 154-023, 
and NGC 784, show indications of distributed SF. These 
findings are confirmed in the maps of the BHeB stellar 
distributions in the Appendix which show that sites of 
recent SF are highly centralized in some galaxies but can 
also extend over much of the optical disk of a galaxy in 
others. The majority of the systems lie somewhere be- 
tween these two extremes. While delineating between 
centralized and distributed SF would be based on arbi- 
trary thresholds, it is clear that not all starbursts are 
highly centralized in the most recent 100 — 250 Myr. 

In Figure |12[ we present a comparison of two con- 
centration metrics calculated over the duration of the 
burst events, which is longer than the 100 — 250 Myr 
timescale in all of the galaxies. Note the plotted range 
for the ratio in birthrate parameters is smaller in Fig- 
There 



12 than in the companion plot in Figure 11 



is a general correlation between the two metrics. When 
the birthrate parameter in the outer regions is within a 
factor of two of the central regions, the fraction of stellar 
mass formed in the central regions is also smaller. One 
galaxy, UGC 6456, is an outlier from the trend. 

Comparing Figures [TT] and |12[ we find a range in the 
spatial extent of the starbursts over both timescales stud- 
ied. This distribution in the spatial extent of the star- 
bursts did not correlate with absolute magnitude, star- 
burst duration, or average SFRs. Note that our sample 
was selected based the HST o bserv ations available in the 
public archive, as discussed in §2.1[ and it is therefore not 
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Figure 11. Comparison of tliree metrics evaluating ttie spatial 
distribution of SF over the most recent ~ 100—250 Myr (see Tablell] 
for time bins on the individual galaxies) . The red "C" in each paneT 
marks the centrally concentrated end of the spectrum, while the 
red "D" marks the spatially distributed end. The three metrics are 
correlated and show a continuous spectrum in the spatial extent of 
SF in these starburst galaxies. 



a volume-limited sample. Thus, we are unable to mea- 
sure the distribution function of the spatial extent of SF 
in starburst dwarf galaxies as a class of objects. However, 
our result does have implications for understanding the 
trigger mechanisms of the starbursts. Whatever phys- 
ical mechanisms initiate these starbursts must account 
for this range in spatial properties. 

The galaxies that show the highest degree of central- 
ized SF (e.g., NGC 1569, NGC 4449, NGC 5253) are 
also the most well-known, and have been identified as 
starbursts according to numerous and diverse observa- 
tional criteria, including h osting a popul ation of Wolf- 
Rayet ( W-R) stars (Schaer er et al. 1999) or super star 
clusters ( jLarsen eF al.jpOQg), having thres hold infrared 
emission line ratios ( ihornle 



ey et al. 1 12000]), or showing 



significant Ha emission (Lee et al. 2009), and high SFR 



su rface d ensities in their central regions (Kennicutt et 
al.||2005|). These systems lend themselves to being clas 



sihed as present-day starbursts because of the observa- 
tional impact made by their centrally concentrated SF, 
particularly when classified with short timescale star for- 
mation tracers (e.g., Ha emission and W-R stars). In 
contrast, the currently more distributed starbursts (e.g., 
ESQ 154-023, UGC 4483, NGC 4068) do not always 
meet this diverse range of observational criteria and are 
thus less readily identified as starbursts. The range in the 
spatial extent of the starbursts highlights a bias towards 
classifying galaxies with centralized SF as starbursts over 
galaxies with more distributed SF. This bias is compli- 
cated by the fact that there is no absolute threshold or 
criterion for classifying a galaxy as a starburst system. 
Some of the most well-known starbursts in the sample 
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Figure 12. Comparison of two metrics evaluating the spatial dis- 
tribution of SF over the duration of each starburst. Similarly to 
Figure [TT] the red "C" in each panel marks the centrally concen- 
trated end of the spectrum, while the red "D" marks the spatially 
distributed end. Note the x-axis has a shorter range than the com- 
panion plot in Figure [TT] There are number of galaxies that lie 
in the lower left quadrant indicating that SF was more distributed 
ove r th ese longer timescales. Comparing these results with Fig- 
^1] we find a range in the spatial extent of the starbursts over 
ETir 



ure 

both timescales studied 



that have the highest degree of centralized SF also have 
the highest SFRs, suggesting a possible correlation be- 
tween the SFR in a starburst and the spatial extent of 
the SF. To explore this possible connection, we ranked 
the galaxies according to the three concentration met- 
rics and then compared the average of these three rank- 
ings to the peak SFR of the starbursts. The peak SFRs 
(Paper H) were normalized by the absolute B-band lu- 
minosity (in units of L©) to elliminate any bias in the 
SFR due to galaxy size. In Figure [TSj we present this 
comparison with the most centralized starbursts ranked 
with smaller numbers and the most distibuted ranked 
with larger numbers. We find no correlation between 
the spatial extent of the SF and the normalized SFR of 
starbursts. 

4. CONCLUSIONS 

We have used optical imaging of resolved stellar pop- 
tilations obtained from the HST data archive to study 
the spatial distribution of SF in fifteen nearby starburst 
dwarf galaxies. We measure the concentration of SF over 
three different timescales using (1) the location of young 
BHeB stars from the V and I band observations, (2) 
newly derived SFHs for central and outer regions of the 
galaxies, and (3) an estimate of the fraction of stellar 
mass created in the central regions. The three metrics 
show that the spatial extent of SF lies on a continuum in 
these galaxies both over the most recent 100 — 250 Myr, 
and over the life of the starburst events. SF is seen to be 
highly concentrated in a some of the sample (e.g., NGC 
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Figure 13. Comparison of the galaxies ranked by their their spa- 
tial concentration of SF with the peak SFR of the starbursts, nor- 
malized by their absolute B-band luminosity in units of L0 . The 
final rankings were determined by averaging the individual rank- 
ings of the three concentration metrics, from the most centrally 
concentrated to the most distributed. The three rankings were in 
general agreement, spanning a range that varied between 1 — 4. 
The peak SFRs were taken from Paper II and normalized to ellim- 
inate bias in the SFR due to the size of a galaxy. Uncertainties in 
the normalized SFRs are smaller than the individual points shown. 
As seen in the Figure, there is no correlation between the spatial 
extent of the starburst and the SFR normalized by luminosity, nor 
was there a correlation with the total SFR. 



4163, NGC 5253, and NGC 6789), distributed in other 
systems (e.g., NGC 4068, ESQ 154-023, NGC 784), and 
somewhere in between for the remaining galaxies. The 
range in spatial extent in the starbursts did not correlate 
with the peak SFR of the starbursts. 



These new findings on the greater spatial extent of 
these starbursts, combined with the recent results mea- 
suring starbursts to last hundreds of Myr (Paper II) , pro- 
vide a new perspective on the bursting mode of SF. In 
contrast to the often reported or assumed smaller, cen- 
tralized spatial scales and short timescales (10 Myr) of 
starbursts, these results show that starbursts can be a 
global phenomenon. In addition, the range in spatial 
extent of the starbursts, coupled with the different ob- 
servational signatures of centralized versus distributed 
SF, highlights a bias towards classifying galaxies with 
centralized SF as starbursts over galaxies with more dis- 
tributed SF. This bias is complicated by the fact that 
there is no absolute threshold or criterion for classifying 
a galaxy as a starburst system. 

While not a volume limited sample of starburst dwarf 
galaxies, these results provide context for understand- 
ing the characteristics of starbursts in low-mass galaxies. 
Any trigger mechanisms for the elevated levels of SF in 
our sample must account for these characteristics. The 
physical mechanism governing the spatial and temporal 
extent of these starbursts was not investigated, but is 
likely based both on the starburst trigger mechanism and 
the equation of state of the gas, including the degree of 
turbulence present in the ISM, in each system. 
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APPENDIX 

MAPS OF THE RECENT SF TRACED BY BHEB STARS 



In Figure 14 we present the CMDs, BHeB luminosity functions, and maps of the sites o f re cent SF for four 
luminosity or age bins for the entire sample. This analysis follows the prescription described in §2.2| In systems with 
higher metallicity (i.e., higher levels of extinction) and photometric crowding, BHeB stars only as old as 100 Myr can 
be unambiguously separated in the CMDs (e.g., NGC 625, NGC 5253, NGC 1569, NGC 4449), whereas BHeB stars 
up to 250 Myr in age could be separated in other systems (e.g., UGC 9128, UGC 4483). The scaling of the asterisks 
was uniquely determined for each system depending on the range of SFRs measured over the time periods studied; 
the size of the asterisks cannot be compared across systems. The preprint version of this work contains low resolution 
color figures; the published version contains high resolution color figures. 



Spatial Distribution of Star Formation 



13 






o 
o 

o 

Ln 



o 

o 
o 
CM 



Maps of BHeB stars of varying ages on smoothed HST images 

Figure 14. Top left panel: CMD of the galaxy with the region containing BHcB stars used in the analysis marked in blue. Top right 
panel: The luminosity function of the BHeB stars identified in the CMD. Middle and lower panels: The location of the BHeB stars of four 
different age ranges overlaid on smoothed V band HST images of the galaxy. The area of the points is scaled by the SFR enabling a direct 
comparison between age bins. The white box drawn on the first image encompasses the area used in reconstructing the SFH of the central 
region. The region outside of this box comprises the outer region used in our analysis. Note that the HST images are not shown in WCS 
orientations. UGC 9f 28 hosts a fossil burst that ended in the last < 15 Myr. 




Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. Note that the WFPC2 PC chip was positioned on the center of the galaxy, thus the analysis was performed within the 
"chevron" shape of the field of view. 




Maps of BHeB stars of varying ages on smoothed HST images 



Figure |14[ Cont. NGC 4163 hosts a fossil starburst that ended ~ 100 Myr ago. Thus, the sites of SF shown in the middle panels (i.e., 
t< 50 Myr and 100<t<50 Myr) are not considered part of the starburst event. We show all four time bins for completeness, but only the 
spatial distribution of the SF between 100 — 200 Myr is considered in our analysis. 
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Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. Note that the WFPC2 PC chip was positioned on the center of the galaxy, thus the analysis was performed within the 
"chevron" shape of the field of view. 




Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. NGC 6789 hosts a fossil starburst that ended ^ 80 Myr ago. Thus, the sites of SF shown in the < 80 Myr panel 
arc not considered part of the starburst event. We show all four time bins for completeness, but only the spatial distribution of the SF 
between 80 — 250 Myr is considered in our analysis. Also note that the WFPC2 PC chip was positioned on the center of the galaxy, thus 
the analysis was performed within the "chevron" shape of the field of view. 




Maps of BHeB stars of varying ages on smoolhed HST images 

Figure |14[ Cont. 




Maps o1 BHeB sftars of varying ages on smoothed HST images 

Figure |14[ Cont. Note the large asterisks reflect a sustained, high level of SF over the time bins studied. 




Maps of BHeB stars of varying ages on snnoothed HST images 

Figure [l4) Cont. 
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Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. Note the large asterisks reflect a sustained, high level of SF over the time bins studied. 
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Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. Note that while NGC 625 hosts a fossil burst than ended < 35 Myr ago when averaged over the entire field of view, 
the central regions in this system show ongoing levels of SF. Thus, we include the SF over the most recent 35 Myr for completeness in the 
analysis of this starburst event. 
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Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. 






Maps of BHeB stars of varying ages on smoothed HST images 

Figure |14[ Cont. 
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Table 1 

The Galaxy Sample and BHeB Star Timescales 







D25 Area 


Deproj. D25 Area 


Duration 


Max. Meas. Age of 


Galaxy 


f mae) 


in FOV (%) 


in FOV (kpc^) 


of Burst (Myr) 


BHeB stars (Myr) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


UGC 9128 


-12.45 


100% 


1.2 


1300±300 


250 


UGC 4483 


-12.68 


100% 


1.0 


>810±190 


250 


NGC 4163 


-13.75 


100% 


3.2 


460±70 


200 


UGC 6456 


-13.85 


100% 


3.0 


>570±60 


200 


NGC 6789 


-14.60 


100% 


1.7 


480±70 


250 


NGC 4068 


-14.96 


100% 


19.5 


>459±50 


200 


DDO 165 


-15.19 


100% 


18. 


>1300±300 


158 


IC 4662 


-15.39 


100% 


2.2 


>450±50 


158 


ESQ 154-023 


-16.21 


63% 


120 


>450±50 


158 


NGC 2366 


-16.33 


96% 


40. 


>450±50 


158 


NGC 625 


-16.26 


63% 


35. 


450±50 


100 


NGC 784 


-16.78 


73% 


64. 


>450±50 


158 


NGC 5253 


-16.98 


100% 


28. 


>450±50 


100 


NGC 1569 


-17.76 


87% 


3.0 


>450±50 


100 


NGC 4449 


-18.02 


87% 


40. 


>450±50 


100 



Note. — Column (1) Galaxy name. Column (2) B— band lumi nosity of the galaxy. Column (3) 
Indicates the percentage of the D25 area ( |Karachentsev et al. | [2004[ | in the observational field of view 
(F.O.V.). Column (4) Measures the deprojected area in the F.O.V. out to the D25 area. Column (5) 
The starburst durations reported in Paper II. Lower limits indicate that starbursts are on-going in these 
systems. Total starburst durations are measured in 4 fossil burst systems. Column (6) The timescale 
over which it is possible to unambiguously select the BHeB from the optical CMD with a maximum age 
of 250 Myr. The maximum age of 250 Myr was chosen as this timescale approximates the lower limit 
for which stellar structure has been shown to exist in dwarf galaxies (Bastian et al. 2011). Diffusion of 
stars from their birth sites becomes more important at longer timescales, thus reducing the accuracy of 
tracing sites of recent SF with the distribution of HeB stars older than ~ 250 Myr. 



Table 2 

The Spatial Extent of the Bursts 



Most Recent ^ 100 - 250 Myr Lifetime of Burst 

Concentration Ratio of Ratio of Central to Outer Fraction M* in Ratio of Central to Outer Fraction M* in 
BHcB/RGB stars (%) Birthrate Parameters Central Region Birthrate Parameters Central Region 



Galaxy 


(2) 


(3) 


(4) 


(5) 


(6) 


UGC 9128 


44 


4.3 ± 0.4 


69 ± 7 


1.3 ± 0.1 


40 ± 4 


UGC 4483 


71 


4.0 ± 0.4 


46 ± 4 


1.6 ± 0.1 


25 ± 5 


NGC 4163 


21 


0.0 ± 0.0 


96 ± 10 


1.3 ± 0.1 


70 ± 13 


UGC 6456 


20 


16.0 ± 2.3 


58 ± 3 


4.8 ± 0.5 


29 ± 6 


NGC 6789 




16.5 ± 0.1 


85 ± 3 


4.3 ± 0.1 


56 ± 13 


NGC 4068 


97 


2.1 ± 0.1 


51 ± 2 


0.1 ± 0.1 


37 ± 1 


DDO 165 


65 


2.5 ± 0.1 


46 ± 2 


1.3 ± 0.1 


32 ± 2 


IC 4662 


25 


5.8 ± 0.1 


78 ± 3 


2.7 ± 0.1 


59 ± 8 


ESO 154-023 


100 


0.9 ± 0.1 


27 ± 2 


1.1 ± 0.1 


30 ± 4 


NGC 2366 


85 


1.5 ± 0.2 


61 ± 6 


1.0 ± 0.1 


52 ± 4 


NGC 625 


10 


4.8 ± 0.1 


63 ± 5 


3.0 ± 0.1 


51 ± 9 


NGC 784 


54 


1.3 ± 0.1 


25 ± 1 


0.8 ± 0.1 


19 ± 2 


NGC 5253 


13 


0.0 ± 0.0 


94 ± 6 


4.0 ± 0.1 


59 ± 6 


NGC 1569 


15 


12.1 ± 0.1 


83 ± 4 


4.7 ± 0.1 


67 ± 8 


NGC 4449 


28 


7.6 ± 0.1 


84 ± 4 


3.4 ± 0.1 


69 ± 13 



Note. — Column (1) Galaxy name. Columns (2) — (4) Results are based on the SF from Myr to 100 - 250 Myr ago (sec Tabic [l] Column 
6) except for NGC 4163 and NGC 6789. In these two post-burst galaxies, the SF was measured from 100 - 200 Myr and 80 - 250 Myr ago 
respectively, corresponding to time periods of actively bursting SF. Column (2) Comparison is based on the number density of the BHeB to 
the RGB stars. Uncertainties are estimated to be of order 10%. Columns (5) — (6) Results are based on the SF that has occured over the 
duration of the starburst in each galaxy. 
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